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Optical tweezers are often applied to control the dynamics of objects by scanning light. However, there is a 
limitation that objects fail to track the scan when the drag exceeds the trapping force. In contrast, 
Laguerre-Gaussian (LG) beams can directly control the torque on objects and provide a typical model for 
nonequilibrium systems such as Brownian motion under external fields. Although stable "mid-water" 
trapping is essential for removing extrinsic hydrodynamic effects in such studies, three-dimensional 
trapping by LG beams has not yet been clearly established. Here we report the three-dimensional off-axis 
trapping of dielectric spheres using high-quality LG beams generated by a special holographic method. The 
trapping position was estimated as ~ half the wavelength behind the beam waist. These results establish the 
scientific groundwork of LG trapping and the technical basis of calibrating optical torque to provide 
powerful tools for studying energy- conversion mechanisms and the nonequilibrium nature of biological 
molecules under torque. 

Optical tweezers have become essential tools for the experimental study of single-molecule biophysics 1 " 3 
and nonequilibrium statistical physics 4 " 7 due to the controllability of the pico -Newton force with nano- 
metre-scale position precision. The basic function of optical tweezers is to trap dielectric objects in three- 
dimensional space where the scattering and gradient forces of light are balanced. This mechanism is often called 
gradient-force trapping 8,9 . It is also possible to control the dynamics of trapped objects by scanning light 4 " 7 ; 
however, this approach is not directly applicable to the object's motion. The trapped object fails to track the scan 
when the viscous drag exceeds the trapping force for extremely quick scans. Another way of controlling dynamics 
is to introduce force fields, e.g. optical vortices (OV) 10 " 18 , which have the special property of carrying orbital 
angular momenta, can produce nearly ideal torque fields producing endless revolutions of dielectric objects 19 " 23 . 
Optical tweezers with OV are expected to provide powerful tools for the experimental studies of nonequilibrium 
statistical mechanics 24 , e.g. nonadiabatic transitions among nonequilibrium stationary states 25 . 

Since the principle of gradient-force trapping is simple and universal, compatibility between optical tweezers 
and OV appears straightforward. This impression, however, might introduce confusion to a general description 
of optical trapping by OV. The term "trapping", which Ashkin introduced to indicate the three-dimensional 
version for Gaussian trapping 8,9 , is now overused for the two-dimensional case in OV trapping 14,15,21,22,26 . So far, 
experimental studies of three-dimensional OV trapping have been reported in the context of equivalence between 
the spin and orbital angular momenta of Laguerre-Gaussian (LG) beams 10 " 13 , which have received general 
attention for the past two decades as a series of the most common and natural embodiments of OV as well as 
the fundamental elements of light propagation 16 . In these experimental reports, the transversal properties of LG 
trapping were studied by observing objects rotating around their own axes at the beam centre, but the axial 
trapping property failed to be firmly established. The properties of LG trapping have also been discussed from a 
theoretical point of view by several groups 21,27 " 32 in relation to the detailed study of light's scattering force. 
Nevertheless, the situation resembles that in the experimental studies except that a few studies discussed the 
axial stability of LG trapping 27,29,32 in on -axis trapping. 

The axial property of LG trapping becomes more obscure in the case of off- axis trapping. It is also experi- 
mentally difficult to achieve three-dimensional off-axis trapping by only a single LG beam of micrometre-sized 
dielectric objects. One reason for the difficulty is that the gradient force in LG trapping tends to become smaller 
than the usual Gaussian trapping due to the lower spatial concentration of the light-power distribution for LG 
beams. Moreover, a common phase-only-modulation holographic method for generating LG beams further 
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Figure 1 | Outlines of experiments and analysis, (a) Optical system. Blue inset shows schematic picture of a submicrometre- sized polystyrene sphere 
trapped by balance between scattering and gradient forces of a LG beam, (b) Data processing flow for deriving a revolution radius from series of sphere 
images. Images of revolving polystyrene spheres were converted to a time series of centroid information. Centre of revolution was first determined by 
fitting a series of centroid data to a circular orbit. Then the radial-position distribution was obtained with respect to the centre, and R s was derived by 
fitting the radial distribution to a Gaussian profile. 



reduces the axial light-power concentration because of the mixture of 
radially higher-order LG-mode contents 33 . Radially higher-order LG 
modes have multi-ringed patterns, where the inner rings are enlarged 
in the axial direction to produce smaller axial gradient forces than the 
lower-order modes. As a result, such additional tricks as gravity 19 , the 
surface of a cover glass 21 ' 22 ' 26 , and another Gaussian trapping beam 34 
are necessary to fix the axial trapping position in the LG trapping. On 
the other hand, the off- axis revolution of small dielectric objects on a 
three-dimensional circular trajectory has recently been achieved 35-37 
by introducing another optical technique. However, this technique 
depends on artificial optical fields specific to the three-dimensional 
revolution, and a fundamental question in optical science remains: 
can a single pure LG beam achieve three-dimensional off-axis trap- 
ping of small dielectric objects? 

This paper reports the stable off- axis trapping of micrometre -sized 
dielectric spheres using holographically generated high-quality LG 
beams under microscopy. A holographic complex- amplitude modu- 
lation scheme 33 was introduced to generate various LG beams of 
mode purities over 95% 38 . This approach enabled three-dimensional 
trapping with sufficient trapping stability for quantitative study. 
Evidence of three-dimensional LG trapping was also established by 
observing the revolutions of micrometre- sized dielectric spheres 
trapped in "mid-water", i.e., completely surrounded by water, not 
touching the bottom or any wall of the chamber. We confirmed 
three-dimensional LG trapping without any LG-independent sup- 
ports due the revolution radius of a dielectric sphere trapped in a 50- 
^m-thick water layer being independent of the LG-beam power. 
Moreover, the axial trapping position was estimated to be ~ half 
wavelength behind the beam waist from the changes in the revolution 
radii; these changes were caused mechanically by pressing the revolv- 
ing spheres onto a pair of glass walls separated by 3 /mi. From 
another viewpoint, the hydrodynamic surface coupling effect of 



viscous drag on the trapped sphere was also simultaneously evalu- 
ated to confirm the trapping in mid-water. 

Since the forces applied to objects are usually derived according to 
the viscous drag of the surrounding medium in microscopy, stable 
three-dimensional LG trapping is essential for determining the abso- 
lute optical torque free from surface effects on viscous drag 39 ' 40 and 
for establishing the technical grounds of optical tweezers as manip- 
ulation tools in the microscopy world. Optical tweezers with LG 
beams are also expected to enhance progress in such fields as the 
torsional properties of biological molecules 41 ' 42 , microrheology 39,40 , 
and nonequilibrium statistical mechanics 4 " 7,43,44 . 

In the following, we restrict our interest to the radially lowest- 
order LG beams that present single-ringed doughnut patterns and 
carry the orbital angular momenta of ifi per photon, where / is an 
integer number referred to as an azimuthal-mode index. A dielectric 
object irradiated by light receives a total force equal to the sum of the 
scattering and gradient forces acting on each small portion of the 
object 45 . Following this simple picture, higher-index dielectric 
objects under the irradiation of an LG beam are expected to be 
revolved on the ring pattern by being trapped on a plane slightly 
behind the beam waist, where the scattering and gradient forces are 
balanced (blue-lined inset in Fig. la). 

Results 

Dependence of revolution radius on LG-beam power. To confirm 
that micrometre- sized dielectric objects can be trapped in mid- water 
on the bright rings of LG beams without any light-independent 
support, we first examined the behaviour of the revolution radius 
(R s ) against incident LG-beam power (P LG ) under the condition 
where a glass surface lies too far from the trapped objects to 
support trapping. Since LG beams keep similar beam patterns of 
different sizes during propagation 16 , the cross-sectional light-power 
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Figure 2 | Behaviour of revolution radius under change of light power for 
polystyrene spheres of various radii, (a), (b), and (c) show experimental 
results for polystyrene spheres of radius 0.2285, 0.3855, and 0.51255 /mi, 
respectively. Errors in revolution radii are evaluated as standard deviations 
of ten measurements but five times smaller than plot symbols. In each 
figure, classifications of trapping LG beams are distinguished by different 
colours and symbols: red, blue, and brown represent / = 2, 3, and 4, 
respectively, and closed and open circles denote left and right 
polarizations. Straight lines show most suitable linear models for the light- 
power dependence of the revolution radius under individual conditions of 
trapping LG beam and sphere radius. Each value denoted as A AIC is a result 
of the subtraction of AIC for a square-root model from that for a linear 
model, providing a benchmark for the goodness of the linear model. 

density should be inversely proportional to the square of the beam 
radius. We start by considering a situation where a dielectric sphere 
was trapped by the balance between the scattering and other light- 
independent forces (e.g., gravity) at P LG . When P LG is increased a 
times, the balance position will shift to where the light-power density 
is about a" 1 times, because the total scattering force is approximately 
proportional to an integral of the light-power density over the 
spherical volume. According to the above mechanism, ,R S must be 

1/2 

proportional to P LG for trapping with light- independent supports, 
while ,R S is independent of P LG for gradient-force trapping. 

Figures la and lb show schematic pictures of the setup and the 
procedure for observing the motions of polystyrene spheres trapped 
in mid-water by holographically generated LG beams 33 (See 
Methods). The changes in R s versus P LG were measured for various 



conditions that were determined as combinations of three different 
sphere sizes, polarizations, and azimuthal-mode index / of the LG 
beams. Here, the values of P LG were measured in advance at the 
sample position under dry conditions. Figures 2a-2c exhibit the 
change of ,R S against P LG , where the R s values agree with the effects 
of polarization and / on the bright- ring radius of LG beams under 
tight-focusing conditions 46 . The revolutions were highly stable and 
uniform around the bright ring of the LG beams, and all plot symbols 
in Figs. 2a-2c cover intervals of over five times the standard deviation 
of measured R s , except the case of trapping by the left-polarized 
/ = 2 LG beam in Fig. 2c, which became unstable and the standard 
deviation became comparable to the plot symbols. 

In contrast, trapping of the 0.5125-/im-radius beads by the left- 
polarized / = 2 LG beam belonged to an on- axis trapping condition, 
where the size of the trapped object exceeded the diameter of the 
trapping LG beam. The stability of on- axis trapping was discussed in 
detail from the viewpoint of deterministic dynamics in Ref. 32 and 
assigned to the competition between the dynamics and the damping. 
In fact, ,R S was not always determined as a well-defined value or it 
became anomalously small compared to that for stable trapping even 
if it were successfully determined, indicating the instability of on -axis 
trapping. On- axis trapping phenomena also suffer from thermal 
fluctuations presenting more complicated dynamics, which is 
beyond the scope of the present study but an interesting topic for 
future study. 

We compared the goodness of the linear and square- root models 
for the P LG -dependence of the ,R S in the context of Akaike's informa- 
tion criterion (AIC) that enables comparisons between two models of 
different degrees of freedom. Although small but meaningful nega- 
tive inclinations were observed in the fitted linear models, the AICs 
for the linear models are at least 36.4 smaller than those for the 
square-root models. Such differences in AIC over 10 X of the model's 
degree of freedom indicate that the behaviour of R s almost comple- 
tely mismatches the assumption of dominant light-independent 
forces. 

Estimation of trapping position. We establish direct evidence of 
LG-beam trapping in mid-water by estimating the axial trapping 
position with an approach that is effective even in regimes smaller 
than the wavelength of light. A focused LG beam, whose behaviour 
resembles a Gaussian beam, shrinks toward the beam waist during 
propagation (we assume the propagation into the z direction) and 
diverges after passing the beam waist position. Using this property, 
the axial trapping position can be estimated by the following idea 
(Figs. 3a-3c). We start from a situation where a polystyrene sphere of 
radius a is trapped in mid-water under the balance between the 
gradient and scattering forces as in Fig. la. In this situation, the 
sphere revolves on the bright ring of the LG beam at trapping 
position z t (dashed line in Fig. 3b). Here, the trapping position 
indicates a relative position in the sample, although z t should 
almost be fixed during the mid-water trapping if measured from 
the beam waist. When the sample is scanned toward the in- 
direction and the sphere reaches the front cover glass (Fig. 3a), ,R S 
will increase while keeping the trapping position z x during the 
"forward" scan. In contrast, when the sample is scanned against 
the z- direction and the sphere reaches the backward cover glass 
(Fig. 3c), the trapping position becomes fixed at z 2 . In this 
situation, ,R S decreases as the trapping position approaches the 
beam waist and starts increasing after passing over the minimum 
radius. Following the above mechanism, we obtain an R s profile if R s 
is measured along with a continuous scan. 

The trapping position was estimated by using a 3 -^m- thick flow 
cell (See Methods). We fixed the polarization of the incident LG 
beams to the right polarization to suppress the relative errors in R s 
and chose P LG as 4.15, 4.38, and 4.00 mW for / = 2, 3, and 4 LG 
beams, respectively, to perform experiments under stable revolution. 
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Figure 3 | Schematic picture of estimating revolution radius by sample scan. Three different situations in sample scan: (a) two-dimensional trapping on 
front interface, (b) three-dimensional trapping free from interfaces, and (c) two-dimensional trapping on back interface. Here, z x indicates the 
position behind the front glass surface by distance a, and z 2 indicates the position in front of the backward glass surface by a. 



Figures 4a-4c show the changes in the measured values according to 
the relative z-position of the sample (in optical length: i.e., scanning 
distance from the minimum R s position multiplied by the refractive 
index of water) for LG beams of / = 2, 3, and 4. Here, red plots show 
the ,R S with horizontal and vertical errors, where the horizontal errors 
are determined as a fifth of the minimum scale of the micrometre for 
sample scanning and the vertical errors are derived in a similar way as 
for the case of Figs. 2a-2c. The R s profiles match the predictions from 
the above "pressing scheme". The plateaus in the R s profile indicate 
trapping in the mid-water, and the plateau widths are commonly 
consistent with \z 2 — Zi\ ( = 2.348 ± 0.166 fim from the sample 
specifications). 

We also evaluated the viscous drag imposed on a revolving sphere 
in each condition to estimate the hydrodynamic coupling effect with 
glass surfaces in terms of Faxen's correction 40 , which appears to 
enhance the viscous drag coefficient around the surfaces (See 
Supplementary Information). The blue plots in Figs. 4a-4c show 
the changes in the viscous drag coefficient that were derived from 
the correlation analyses on the azimuthal dynamics of the dielectric 
spheres. In fact, the small viscous drag domains are consistent with 
the plateaus in the R s profile, which also supports the trapping in 
mid-water between the glass surfaces. Following the pressing 
scheme, the trapping position roughly corresponds to the "back- 
ward" edge of the plateau, which is commonly about 0.6 times the 
light wavelength behind the beam waist for / = 2, 3, and 4. Although 
this is not a quantitative estimation for the trapping position, it is 
nevertheless the first direct evidence for three-dimensional off-axis 
trapping with LG beams. 

Discussion 

In the present study, we estimated the trapping position by scanning 
the samples in the propagation direction of the trapping LG beams. 
Although the scan range was at most a few wavelengths, this can 
change the incident optical condition of the LG beams to the glass/ 
water interface. For example, subtle inclinations are observed at the 
plateaus of R s in Figs. 4a-4c, which can arise from the scanning effect. 
Nevertheless, we recognise it as minor in the present results because 
the LG-beam shapes around the beam waists are reproduced in 
Figs. 4a-4c. This means that the present scanning scheme is also 
effective for roughly observing the propagation behaviour of the 
LG beams under tightly-focused conditions. Even though the heat 
transfer from the trapping LG beams to the sample is also supposed 
to affect the refractive index of water 47 , the change is estimated to be 

0.01%/°C, which will produce only minor optical effects. 

On the other hand, the viscous drag coefficients in Figs. 4a-4c 
present odd behaviour in the z > 0 domain, although the viscosity 
is expected to increase monotonously here. The main origin of this 
phenomenon might be the LG beams reflected from the backward 
glass surface. In the z > 0 domain, the dielectric spheres usually 
revolve around the backward surface (Fig. 3c) and suffer from the 
effects of the reflected LG beams. The scattering force of the reflected 
light reduces the total force that presses the spheres onto the glass 
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Figure 4 | Revolution radius and viscosity versus sample position. 

Behaviour of revolution radius (red plots: labeled on left axis) and viscous 
drag coefficient (blue plots: labeled on right axis) are shown as functions of 
sample position for (a) / = 2, (b) / = 3, and (c) / = 4 LG beams. Here, 
vertical error bars indicate standard deviations of ten measured values. 
Each sample position in horizontal axis is as sample-scan distance 
measured from "beam waist" position (sample position marking 
minimum revolution radius) in optical length, i.e., relative scan distance 
multiplied by the refractive index of water. The error in sample position is 
determined as a fifth of the minimum scale of micrometre for sample 
scanning. Hatched areas indicate mid-water trapping regimes 
corresponding to situation of Fig. 3b to guide the eyes for comparing 
behaviours of revolution radius and viscous drag. 
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surface, and they can move away from the surface to present a super- 
ficial decrease of Faxen's correction on the viscous drag coefficient. 
This hypothesis also matches the absence of similar behaviour in 
z < —2.52, where the objects are likely to revolve around the front 
glass surface and the effects of the reflected LG beams are ignored. 

The present setup, which essentially combines a high-speed 
CMOS image sensor and a computer- controllable spatial light 
modulator (SLM), is expected to provide a powerful tool for the 
experimental study of nonequilibrium statistical mechanics. In par- 
ticular, the study of the information-to-heat conversion 43,44 , which 
requires the detection of an object's position and the switching of a 
force field based on the observed position information in the milli- 
second regime, is further developed by the present experimental 
technology, if the response time of SLM devices is further improved. 

Methods 

Samples. We suspended three sizes of polystyrene spheres (Polysciences Inc., 
Polybead® Cat# 07310, Cat# 07309, and Cat# 07307, diameters 1.025 ± 0.010, 0.771 
± 0.024, and 0.457 ± 0.011 fim, respectively) in ultra pure water (MilliQ) and put 
between 24 X 40 and 18X18 mm 2 cover glasses (NEO, Matsunami: No. 1) with 50- 
/im-thick spacers. The edges of the smaller cover glasses were sealed with resin to 
make flow cells for the measurements of R s versus various P LG . For the trapping- 
position estimation, Cat# 07309 polystyrene spheres were infused with a few 3-fim- 
diameter spheres (Polybead® Cat# 17134: diameter 3.119 ± 0.164 fim) that will 
eventually become spacers. 

Optical setup. LG beams were generated from linearly polarized laser light 
(Millennia® IR, Spectra Physics; wavelength X = 1,064 nm) as first-order diffraction 
components from holograms displayed on an SLM (XI 0468, Hamamatsu Photonics). 
The generated LG beams were introduced to an inverted microscope (1X71, 
Olympus) equipped with a stable mechanical stage (KS-O, Chuukousya Seisakujo) 
and focused on samples through an oil-immersion objective (UPLFLN 100X UIS2, 
numerical aperture 1.30, Olympus) after being converted to circular polarization with 
a quarter-wave plate. The visible images of the polystyrene spheres were observed by a 
high-speed CMOS image sensor (C8201, Hamamatsu Photonics), which can acquire 
and process 6,337 images per second at the cost of a limited active area (64 X 64 pixels, 
each of which covers a 0.1 X 0.1 fim 2 area on the observing plane). The image sensor 
was mounted on a position-adjustment rail to observe the images under off-focus 
conditions. Trapping positions were estimated by scanning the sample stage in one 
way of the z direction, manually turning a micrometre and adjusting the sensor 
position. The change in the image magnification due to the sensor-position 
adjustment was estimated to be less than ±0.25% with a test pattern ruled into 10-/mi 
squares, meaning that errors in R s are less than statistical deviations. Notice that 
precise adjustment and selection of the optical elements is essential for trapping with 
LG beams. 

Hologram design. We applied a holographic complex- amplitude modulation 
method to generate high-quality LG beams 33 ' 38 of arbitrary orders. The holograms 
were designed by the Kirk- Jones method, and four pixels were chosen for the blazing 
pitch. The beam- size parameters were determined to obtain LG beams that have a 
common Gaussian radius of 1.42 mm at the objective's entrance pupil. The 
holograms were displayed on the SLM after adding a compensation pattern for the 
optical distortion and a Fresnel lens pattern (radius of curvature = 2,759 mm) for 
adjusting the focus position of the X = 1,064 nm LG beams. 

Analysis. Images of the revolving polystyrene spheres were converted to a time series 
of centroid information. Typically, 40,000 images were captured in a single 
measurement, and the measurements were repeated ten times under the same 
conditions to compose a sequence of experiments. The centroids were determined 
with a position uncertainty less than 0.01 pixels corresponding to 1 nm under the 
present setup. The revolution's centre was first determined by fitting a series of 
centroid data to a circular orbit. Then the radial-position distribution was obtained 
with respect to the centre, and the R s was derived by fitting the radial distribution to a 
Gaussian profile (Fig. lb). Since the fitting errors in R s were typically 0.1 nm under 
mid-water trapping conditions, the total errors in R s were dominated by statistical 
deviation. 
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